A new determination of the Pomeron intercept in hard processes 
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A method allowing for a direct comparison of data with theoretical predictions is proposed 
for forward jet production at HERA. It avoids the reconstruction of multi-parton contributions 
by expressing the experimental cuts directly as correction factors on the QCD forward jet cross- 
section. An application to the determination of the effective Pomeron intercept in the BFKL-LO 
parametrization from da/dx data at HERA leads to a good fit with a significantly higher effective 
intercept, otp = 1.43 ± 0.025(stat.) ± 0.025(syst.), than for proton (total and diffractive) structure 
functions. It is however less than the value of the pomeron intercept using dijets with large rapidity 
intervals obtained at Tevatron. We also evaluate the rapidity veto contribution to the higher order 
BFKL corrections. The method can be extended to other theoretical inputs. 

' 1- Introduction 

o : 

The study of forward jets at colliders is considered as the milestone of QCD studies at high energies, since it 
provides a direct way of testing the perturbative resummations of soft gluon radiation. More precisely, the study of 
one forward jet (w.r.t. the proton) in an electron-proton collider 0] seems to be a good candidate to test the energy 
dependence of hard QCD cross-sections. It is similar to the previous proposal of studying two jets separated by a large 
1 rapidity interval in hadronic colliders [Q] , for which only preliminary results are available || . This test is also possible 
p^j | in 7*-7* scattering jjj but here the statistics and the energy range are still insufficient to get a reliable determination 
of the physical parameters for hard QCD cross-sections. Indeed, the proposed (and favored for the moment being) 
set -up jlj is to consider jets with transverse momentum kx of the order of the photon virtuality Q allowing to damp 
the QCD evolution as a function of kx (DGLAP evolution Q) in favor of the evolution in energy at fixed kx (BFKL 
evolution ||). 

Since proposal [jjj was made, a set of interesting studies have been performed to check its relevance. On the 
^v^j ■ experimental ground, HI |Q and ZEUS Js) have published useful results with appropriate cuts (to be displayed later 
on) at relatively small x. On the theoretical ground, the general formulation and some quantitative estimates have 
been performed prior to experiments Q confirming the interest in such processes. The recent theoretical analyses 
have been mainly based on the use of Monte-Carlo simulations, including the multi-parton cross-sections and starting 
from the various frameworks in competition |Io(| . Quite a few analyses arrive at a satisfactory description of the data, 
O { taking into account the specific parametrizations which are choosen. Indeed, the BFKL-based Monte-Carlo (K ] lead 
to quite satisfactory results, while those based on DGLAP evolution meet some difficulty to describe the dataF. 

However, there still remains a problem in the interpretation of those results. Due to the difficulty in handling the 
experimental cuts without introducing in the simulation the whole set of theoretical n-parton contributions to the 
cross-section, it appears difficult to avoid the uncertainties of the reconstruction (with the parameters and constraints 
which are needed to define the scheme in practice). It seems thus difficult to determine unambiguously genuine 
theoretical parameters defining the cross-section one is looking at. One example is the dependence in parameters such 
as infra-red cut-offs, which are not a-priori required in the expression of the total da/dx jet cross sections. Another 
illustration is the so-called "consistency constraint" which appears very useful |ll|] in the expression of next-to-leading 
corrections to the BFKL formula coming from the n-parton contributions, but again is not expressed in terms of the 
da/dx jet cross sections itself. In fact, it does not seem easy, in those schemes, to extract with some precision the 
value of the effective Pomeron intercept ap, i.e. the main theoretical parameter describing the theoretical energy 
dependence in this process. As we know, this parameter is of primordial importance to evaluate the amount of 
next-leading corrections in a BFKL framework jl2) and to confront its effective value with the recent theoretical 
determinations Jl3| . 

We want to address this problem in a quite different way, that is on focussing on the jet cross section da/dx ob- 
servable itself, by a consistent treatment of the experimental cuts and minimizing the uncertainties for that particular 
observable. Let us remark that our approach is not intended to provide a substitution to the other methods, since 
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1 Note however, that some more refined versions of DGLAP evolution including contributions from the resolved off- mass- shell 
photon can describe the data JlOJ|. 
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the Monte-Carlo simulations have the great merit of making a set of predictions for various observables. Hence, our 
method has to be considered as complementary to the others and dedicated to a better determination of the effective 
Pomeron intercept using the da/dx data. As we shall see, it will fix more precisely this parameter, but it will leave 
less constrained other interesting parameters, such as the cross-section normalization. 

One fruitful outcome of the method proposed in the present paper is the possibility of comparing the effective 
intercept with its determination in other processes involving QCD at high energy. In fact, using the parameters 
determined from forward jets at HERA, it is possible to compare with double jet production at Tevatron following 
Ref. for which preliminary experimental analyses have been performed ||] and find a high value of the intercept 
(ap = 1.7±.l±.l in Ref. Q). It can also be confronted with the effective BFKL analysis of proton structure functions 
at small- 2:3,, which give rather low values^](ap ~ 1.1 — 1.21 see Ref. JlJ]). However, in those cases, the result may be 
different, since non-perturbative effects related to the "soft" proton scales are expected to influence the determination 
of parameters. 

Thus, the comparison of the effective BFKL parameter ap obtained for the forward jet production cross-section 
allows for a study of QCD at high energy, aiming at a better understanding of the corrections to the leading-order 
BFKL predictions jjjjll]]. 

The plan of our study is the following: in section 2, we introduce the QCD formalism and our method for determining 
ap. In the following section 3 we determine the kinematic correction factors to the forward jet cross-section data on 
da/dx due to the experimental cuts. In the following section 4, we perform and discuss a (separately and then 
common) fit to HI and ZEUS data. This determines the BFKL parameter ap which is subsequently used in section 4 
for a comparison with the two-jet cross-section at Tevatron from the experimental DO analysis. Discussions on these 
results and comparison with the BFKL study of (total and diffractive) structure functions are presented in section 5 
and conclusions and outlook in section 6. 

2. Formalism 

The cross-section for forward jet production at HERA in the dipole model reads H: 



d^a 



TrNca 2 as(k^) 



exp{e( 7 ,0)y} 



dxdQ 2 dxjdk^d^ 

h T {l) + h L ("i) 
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where 



Y = In 



xj 



e(7,p) = a [2ip(l) -ip(p+ 1 -7) - ^0 + 7)] 
feff(x,n}) = G(x^ 2 ) + ~E(Q, + Qf) 



(1) 

(2) 
(3) 
(4) 
(5) 



are, respectively, Y the rapidity interval between the photon probe and the jet, e( 7 ,p) the BFKL kernel eigenvalues, 
f e ff the effective structure function combination, and /if the corresponding factorization scale. The main BFKL 
parameter is a, which is the (fixed) value of the effective strong coupling constant in LO-BFKL formulae. Note 
that we gave for completion the full BFKL formula including the azimuthal dependence but we will stick to the 
azimuth-independent contribution with the dominant exp{e(7, Q)Y} factor. 
The so-called "impact factors" 



fix 



a s (k 2 ) (r(i- 7 )r(i + 7)) 3 

3tt7 T(2 - 2 7 )r(2 + 2 7 ) 1 



(l + 7)(l"i) 
7(1 - 7) 



(6) 



are obtained from the Ut factorization properties [fP7f of the coupling of the BFKL amplitudes to external hard probes. 
The same factors can be related to the photon wave functions (R]|l(| within the equivalent context of the QCD dipole 
model fnj. 



Note however that taking into account the full BFKL formula may lead to higher dp, namely 1.2 
total structure functions and even reach .4, see Ref. pa], for diffractive proton structure functions. 



1.3, see Ref. M, for 
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Our goal is to compare as directly as possible the theoretical parametrization (|lj) to the data which are collected in 
experiments [@J§. The crucial point is how to take into account the experimentally defined kinematic cuts listed in 
Table I for the reported three sets of data (two for HI with kx = 3.5 or 5 GeV, and one for ZEUS). 



HI cuts 


ZEUS cuts 


E e > 11 GeV 
160 <6' e < 173 deg. 

y >o.i 

7< dj e t <20 deg. 
krjet > 3.5 or 5 GeV 
x jet >0.035 

0.5 < ^ < 2 
10" 4 < x < 4. 10~ 3 


E' e > 10 GeV 

y>0.1 
9 jet >8.5 deg. 
k Tjet > 5 GeV 
x jet >0.036 

0.5 < g& < 2 
4.5 10" 4 < x < 4.5 10" 2 



Table I- Experimental cuts (Hl/ZEUS) 



The main problem to solve is to investigate the effect of these cuts on the determination of the integration variables 
leading to a prediction for da/dx from the given theoretical formula for d^a as given in formula (jl|). The effect 
is expected to appear as bin-per-bin correction factors to be multiplied to the theoretical cross-sections for average 
values of the kinematic variables for a given x-bin before comparing to data (e.g. fitting the cross-sections). 

The idea of our method is threefold: i) for each x-bin, determining the average values of x, Q 2 , Ej, kj* from a 
known and reliable Monte-Carlo simulation of the cross-sections. For this sake, we use the Ariadne Monte-Carlo 
programme plf ; ii) choosing a set of integration variables over d^ a in (hi) in such a way to match closely the 
experimental cuts and minimize the variation of the cross-sections over the bin size; iii) fixing the correction factors 
due to the experimental cuts for each a;-bin, by a random simulation of the kinematic constraints with no dynamical 
input. 

The point i) proposed already in p0|| allows a determination of which average values of the kinematic variables 
have to be taken in the theoretical formula ([!]) for each experimental s-bin. The point ii) comes from the crucial 
requirement to minimize the variation (over the x-bin) of the variables to be retained for the integration. Indeed, 
since the integration procedure multiplies the central value of the integrand by the size of the integration bins, it is 
compulsory to choose adequate variables which lead to a smooth dependence of the integrand and of the effect of the 
kinematic cuts. 

This double stringent requirement can be solved for the forward jet da/dx. For this sake we choose 

The property of this non-trivial choice is the following. The integration variables are choosen in such a way that the 
expression in the square brackets [Q 6 ...] in (|j) is dependent on the ratio and not on each scale separately. Looking 
at the experimental cuts (see Table I), it becomes clear that the choice of this scale-invariant integrand minimizes the 
variation of the observable on the bin, while each scale k T and Q 2 presents large variations and thus would generate 
large integration errors. Indeed, various numerical studies we have performed have demonstrated that it was a sine 
qua non stability condition for the fits. The overall azimuthal integration (A$ = 2tt) cancels the second term in ([!]). 

3. Correction Factors 

The experimental correction factors have been determined using a toy Monte-Carlo designed as follows. We gen- 
erate flat distributions in the variables k^,/Q 2 , 1/Q 2 , xj, using reference intervals which include the whole of the 
experimental phase-space (the $ variable is not used in the generation since all the cross-section measurements are cj) 
independent). In practice, we get the correction factors by counting the numbers of events which fulfill the experimen- 
tal cuts given in Table I for each x-bin. The correction factor is obtained by the ratio to the number of events which 
pass the experimental cuts and the kinematic constraints, and the number of events which fullfil only the kinematic 
constraints, i.e. the so-called reference bin. 



da 
dx 



dxdQ 2 dxj dki, e?$ 
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x 


(J 


O' 1 


Ej e t 




Corr. Factor (.10 -3 ) 


U.UUUoD 


zUz.D 


lo.y 


on 
OZ.O 


4.0 


U.Z7U 


0.00073 


342. 


21.5 


34.4 


5.0 


0.993 


0.0012 


224. 


26.9 


36.9 


5.5 


1.14 


0.0017 


138. 


31.4 


38.1 


5.8 


1.11 


0.0024 


67. 


38.1 


38.8 


6.3 


0.921 


0.0035 


32. 


47.0 


37.9 


6.9 


0.711 



Table Ha- Average values of kinematic quantities and correction factors - HI kx > 3.5GeU 



X 


a 


Q 


Ej e t 


kx 


Corr. Factor (.10 -3 ) 


0.00036 


27.5 


18.0 


35.9 


5.5 


0.108 


0.00073 


126. 


27.0 


36.5 


5.8 


0.695 * 


0.0012 


132. 


32.2 


37.9 


6.3 


0.895 


0.0017 


96. 


34.8 


39.3 


6.5 


0.979 


0.0024 


55. 


40.1 


39.4 


6.7 


0.870 


0.0035 


28. 


48.2 


39.6 


7.2 


0.696 



Table lib- Average values of kinematic quantities and correction factors - HI kx > 5GeV 



X 


a 


Q 2 


E jet 


kx 


Corr. Factor (.10~ 3 ) 


0.0006 


114.0 


28.0 


36.5 


6.3 


0.304 * 


0.0011 


96.2 


39.0 


38.0 


6.9 


0.656 


0.0019 


77.8 


50.7 


39.9 


7.6 


0.966 


0.0033 


34.4 


75.6 


43.8 


8.7 


0.996 


0.006 


14.1 


113.6 


49.6 


10.4 


0.995 


0.01 


6.53 


176.4 


58.5 


12.9 


0.896 * 


0.018 


2.65 


244.7 


67.3 


15.1 


0.653 * 


0.031 


0.65 


366.8 


78.8 


18.8 


0.373 * 



Table He- Average values of kinematic quantities and correction factors - ZEUS kx > 5GeV 

The correction factors are given in Table Ila for HI (kx > 3.5 GeV), Table lib for HI (kx > 5 GeV), and Table lie 
for ZEUS bins together with the value of the bin centers determined with the full Monte-Carlo simulation , and 
the experimental values of the cross-sections^]. We note that the correction factors are quite different from one x-bin 
to an other and much less than one (in 10~ 3 units), cxplicitely showing that the experimental cuts play an important 
role in the cross-section measurement, and that these factors are compulsory to be taken into account if we want to 
get a direct comparison with the theoretical cross-sections. We also note that the correction factors are very much 
different from one another at very low x, showing that the acceptance of these bins is quite low. This is also why it 
is not so easy to be able to get a correct value of the measured cross-section after cuts in those bins. We also get the 
same order of magnitude for the correction factors for the HI and ZEUS experiments because the experimental cuts 



3 Note that we did not use the full Monte-Carlo to get the correction factors in order to avoid any strong model dependence 
as these factors are only due to kinematic effects. It is however more difficult to use a toy Monte Carlo to get accurate values 
for the bin centers, and this is why we used a full Monte Carlo for this sake. However, the dependence of the theoretical 
cross-section on the bin centers is minimized by our specific choice of kinematic variables (see formula (7)). 
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are quite similar. The differences between both experiments are due mainly to the fact that the range in x and Q 2 is 
much lower for HI than for ZEUS (the reference bin for HI goes to lower Q 2 compared to ZEUS). 

4. Fits 

Using the kinematic correction factors determined as described in the previous section, we perform a fit to the HI 
and ZEUS data with only two free parameters, these are the effective strong coupling constant in LO BFKL formulae 
a corresponding to the effective Lipatov intercept ap = 1 + 41og2aAc/7r, and the cross-section normalisation. The 
obtained values of the parameters and the x 2 °f the fit are given in Table III for a fit to the HI and ZEUS data 
separately, and then to the HI + ZEUS data together. Note that one HI point at kx > 5 GeV (7.3 1CU 4 ), and four 
ZEUS points (x = 4. 10 -4 , and the three highest-x points), were not taken into account in the fit and are distinguished 
in Tables II with a star. We will discuss this selection in a little while. 



fit 


a 


ap 


Norm. 


X Z (/dof) 


HI 

ZEUS 
Hl+ZEUS 


0.17 ± 0.02 ± 0.01 
0.20 ± 0.02 ± 0.01 
0.16 ± 0.01 ± 0.01 


1.44 ± 0.05 ± 0.025 
1.52 ± 0.05 ± 0.025 
1.43 ± 0.025 ± 0.025 


29.4 ± 4.8 ± 5.2 
26.4 ± 3.9 ± 4.7 
30.7 ± 2.9 ± 3.5 


5.7 (/9) 
2.0 (/2) 
12.0 (/13) 



Table III- Fit results 

The x 2 °f the fits have been calculated using statistical error only and are very satisfactory (about 0.6 per point for 
HI data, and 1. per point for ZEUS data). We give both statistical and systematic errors for the fit parameters. The 
values of the Lipatov intercept are close to one another and compatible within errors for the HI and ZEUS sets of 
data, and indicate a preferable medium value (ap = 1.4 — 1.5). We also notice that the ZEUS data have the tendency 
to favour a higher exponent, but the number of data points used in the fit is much smaller than for HI, and the HI 
data are also at lower x. The normalisation is also compatible between ZEUS and HI. The fit results are shown in 
Figure 1 and compared with the HI and ZEUS measurements. 

Let us discuss our selection criterium for the fits. Both lowest x points for HI and ZEUS show large correction 
factors but only the lowest x point for ZEUS lies a bit above the prediction, which shows the relevance of the correction 
factors we determined. On the other hand, the three highest x points for ZEUS cannot be described by a BFKL fit 
probably because the x- value is too high (x > 10~ 2 ). Consider now the second lower x point at kx > 5 GeV for the 
HI experiment that we suppressed from the fit (see Table lib). If we include it in the fit the % 2 value goes from 5.7 
to 32, which is due to the small statistical error of this data point (the systematic error is on the contrary very large). 
By comparison, including the lowest x point for ZEUS changes the x 2 from 2.0 to 7.9. In the same way, including 
the highest x points still increases the x 2 to 67.4, showing clearly that these highest x points cannot be described 
using the BFKL formalism. It is interesting to note that all similar discrepancies appear also in other types of fitting 
procedures, e.g. in Ref. [p2[ . 

5. Comparison with other processes 

The final result of our new determination of the effective pomeron intercept is ap = 1.43 ± 0.025 (stat.) ±0.025 
(syst.). This high value of the intercept leads to the following remarks. Our analysis confirms the trend observed 
using DGLAP based Monte-Carlo jll| which have difficulties to reproduce the forward jet cross-section due to a low 
effective pomeron intercept when both k T and Q 2 scales are of the same order. 

On the other hand, our method allows a direct comparison of the intercept values with those obtained in other 
experimental processes, i.e. 7*7* cross-sections at LEP ELjct-jet cross-sections at Tevatron at large rapidity intervals 
F2 and F2 proton structure function measurements |pT[-^6[ . Let us first consider the known determinations of the 
effective intercepts in F2 and measurements at HERA ]23[ |. It is known that the effective intercept determined in 
these measurements is rather lowQ(1.2-1.3). This is the reason why these data can be both described by a DGLAP or 
a BFKL-LO fit Q. 



4 It is interesting to note that the "hard" Pomeron intercept obtained within the framework of two-Pomeron models [^4| fits 
with our determination. However our parametrization ( |l|) c orresponds to only one Pomeron. 

5 Note that in the BFKL descriptions of these data ]15| , fl6| , the effective intercept is taken to be constant, while the Q 2 
dependence comes horn the BFKL integration (see for instance formula ([j])) 
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FIG. 1. The HI data (k T > 3.5 GeV, k T > 5 GeV), and the ZEUS data are compared with the result of the fit. ZEUS data 
are also displayed in logarithmic scales in vertical coordinates to show the discrepancy at high x values. 

Now let us consider processes initiated by two hard probes which allow a more direct comparison between experi- 
ments and BFKL predictions. These processes suppress DGLAP evolution by selecting; events with comparable hard 
scales for both hard probes. Recent data on 7*7* cross-section measurements at LEP |2j| lead to a BFKL description 
with a low effective intercept compatible with the one of F2 and F® at HERA (ap=1.2-1.3 0) []. The fact that 
similar values of the intercepts are found could be interpreted by sizeable higher order corrections to BFKL equation. 
On the other hand, it is interesting to note that our result based on forward jet measurement at HERA obtained 
in comparable Q 2 (Q 2 ~ 10 GeV 2 ) and rapidity (Y ~ 3-4) domains is quite different. The value of the intercept is 
significantly higher. 

It is also fruitful to compare our results with the effective intercept we obtain from recent preliminary dijet data 
obtained by the DO Collaboration at Tevatron || . The measurement consists in the ratio R = crisoo / &630 where a is the 
dijet cross-section at large rapidity interval Y ~ A77 for two center-of-mass energies (630 and 1800 GeV), Ar7 180 o = 4.6, 
A?763o = 2.4. The experimental measurement is R = 2.9 ±0.3 (stat.) ±0.3 (syst.). Using the Mueller-Navelet formula 
H , this measurement allows us to get a value of the effective intercept for this process F] 



The statistics for these data is still very low. L3 and OPAL Collaborations have released the cuts used to enhance BFKL 
effects to get more statistics [^Q. These data can be both described by BFKL and DGLAP evolution equations. 

7 Formula (^) is obtained after integration over the jet tranverse energies at 630 and 1800 GeV, Et x , Et 2 ■ We note that a non 
integrated formula shows a sizeable dependence on Et x / Et 2 , which could be confronted with experiment [E6[ . 
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We get ap=1.65 ± 0.05 (stat.) ± 0.05 (syst.), in agreement with the value obtained by DO using a saddle-point 
approximation |J. This intercept is higher than the one obtained in the forward jet study. 

The question arises to interpret the different values of the effective intercept. It could reasonably come from the 
differences in higher order QCD corrections for the BFKL kernel and/or in the impact factors depending on the initial 
probes (7* vs. jets). In order to evaluate the approximate size of the higher order BFKL corrections, we will use their 
description in terms of rapidity veto effects p7| . In formula ([{]), we make the following replacement 

exp (e(7 , 0)Y) - E~ 0(Y - (n + 1)6) &h °) ^ + « . (9) 

The Heaviside function 6 ensures that a BFKL ladder of n gluons occupies (n + l)b rapidity interval where 
b parametrises the strength of NLO BFKL corrections. The value of the leading order intercept is fixed to 
a p = 1.75(as(Q 2 = 10) = 0.28), where Q 2 — 10 GeV 2 is inside the average range of Q 2 in the forward jet mea- 
surement. The fitted value of the b parameter obtained using the forward jet data is found to be 1.28 ± 0.08 (stat.) 
± 0.02 (syst.). Imposing the same value of ap with Tevatron data gives 6=0.21 ± 0.11 (stat.) ± 0.11 (syst.). Note 
that the theoretical value of b for the NLO BFKL kernel is expected to be of the order 2.4, which is also compatible 
with the result obtained for the 7*7* cross-section. A contribution from the NLO impact factors is not yet known, 
and could perhaps explain the different values of b. 

6. Conclusions 

To summarize our results, using a new method to disantangle the effects of the kinematic cuts from the genuine 
dynamical values of the forward jet cross-sections at HERA, we find that the effective pomeron intercept is ap — 
1.43±0.025 (stat.) ±0.025 (syst.). It is much higher than the soft pomeron intercept, and, among those determined in 
hard processes, it is intermediate between 7*7* interactions at LEP and dijet productions with large rapidity intervals 
at Tevatron. 

Looking for an interpretation of our results in terms of higher order BFKL corrections expressed by rapidity gap 
vetoes b between emitted gluons, we find a value of b =1.3, which is sizeable but less than the theoretically predicted 
p2| value for the NLO BFKL kernel (b =2.4). The observed dependence in the process deserves further more precise 
studies 

Last but not least, the derivation of the correction factors given in Table II is independent of the theoretical input 
and could be used to test any model suitable for the jet cross-section. 
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